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ABSTRACT 

Context. A substantial fraction of the stellar mass growth across cosmic time occurred within dust-enshrouded environ- 
ments. So far the identification of complete samples of distant star-forming galaxies from the short wavelength range 
has thus been strongly biased by the effect of dust extinction. Yet, the exact amount of star- forming activity that took 
place in high-redshift dusty galaxies currently missed by optical surveys has been barely explored. 
Aims. Our goal is to determine the fraction of luminous star- forming galaxies at 1.5 < z < 3 potentially missed by the 
traditional color selection techniques because of dust extinction. We also aim at quantifying the contribution of these 
sources to the IR luminosity and cosmic star formation density at high redshift. 

Methods. We base our work on a sample of 24/im sources brighter than 80/iJy and taken from the Spitzer survey of 
the COSMOS field. Almost all of these sources have accurate photometric redshifts. We apply to this mid-IR selected 
sample the BzK and BM/BX criteria as well as the selections of the ^^IRAC peakers^^ and the Optically -Faint IR-hrighV^ 
galaxies, and we analyze the fraction of sources identified with these techniques. We also computed 8/im rest-frame 
luminosity from the 24/im fluxes of our sources, and considering the relationships between Ls^m and Lpa^ and between 
Ls^rn and Ljr we derived pm and then psFR for our MIPS sources. 

Results. The BzK criterion offers an almost complete (^^90%) identification of the 24/im sources at 1.4 < ^ < 2.5. On the 
contrary the BM/BX criterion miss 50% of the MIPS sources. We attribute this bias to the effect of extinction which 
redden the typical colors of galaxies. The contribution of these two selections to the IR luminosity density produced 
by all the sources brighter than 80/iJy are from the same order. Moreover the criterion based on the presence of a 
stellar bump in their spectra (''IRAC peakers" ) miss up to 40% of the IR luminosity density while only 25% of the 
IR luminosity density at z^2 is produced by Optically- Faint IR-hrighV^ galaxies characterized by extreme mid-IR to 
optical flux ratios. 

Conclusions. Color selections of distant star-forming galaxies must be used with lots of care given the substantial bias 
they can suffer. In particular, the effect of dust extinction strongly impacts the completeness of identifications at the 
bright end of the bolometric luminosity function, which implies large and uncertain extrapolations to account for the 
contribution of dusty galaxies missed by these selections. In the context of forthcoming facilities that will operate at 
long wavelengths (e.g., JWST^ ALMA, SAFARI, EVLA, SKA), this emphasizes the importance of minimizing the 
extinction biases when probing the activity of star formation in the early Universe. 

Key words. Galaxies: high-redshift - Infrared: galaxies - Cosmology: observations 



1. Introduction 

It is well established that the peak of star-forming activ- 
ity in the Universe and the bulk of stellar mass assembly 
in galaxies occurred at 1 < z < 3 (e.g., Madau et al [ 
T9961 ISteidel et al.[ [19991 IDickinson et al.[ [2003; Ho pkins 



& Beacom 2006 Arnouts et al. , 2007). However the role 



and the contribution of the different processes that gov- 
erned this build-up of stellar mass are still open questions. 
Improving our physical understanding of the population 
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of star-forming galaxies that contributed to the growth of 
structures in the distant Universe thus remains a critical 
issue for modern extragalactic astrophysics. 

To reach this goal, the broad variety of physical prop- 
erties observed at high redshift usually implies the use of 
large and complete samples of galaxies mainly selected as 
a function of star formation rate ( SFR) or stellar mass 



see for instance Elbaz et al. , 2007; Noeske et al. , 2007 



Ilbert et al. 2Q1Q[ ). This requires systematic identincations 
of sourceswith spectroscopic or photometric redshifts (e.g., 
Wolf et all [20031 IFontana et all [2004| jClazebrook et al.l 
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2004 Cimatti et al. 2008). However, the determination 
of reliable spectroscopic redshifts for such large samples 
of distant galaxies is time-consuming because of the dif- 
ficulty to identify emission or absorption lines with high 
enough signal to noise. Flux limited optical selections are 
also dominated by populations of low-to-intermediate red- 
shift sources and they result in statistically small numbers 
of distant objects, while the limited coverage of spectro- 
scopic surveys makes them sensitive to the cosmic variance 
effect. Photometric redshifts, on the other hand, can suffer 
various systematics and contamination from catastrophic 
failures depending on the wavelength range and the num- 
ber of bands they are based on, and their uncertainties can 
become significant in the case of distant galaxies. 

To minimize these difficulties, various techniques based 
on single- or two-color criteria have been proposed to iden- 
tify specific populations of galaxies in the early Universe. 
Among them the first and best-known example is proba- 
bly the selection of Lyman Break Gal axies from the typical 



UnGR colors of sources at z ~ 3 (Steidel et al., 2003) 
Similarly, other techniques have been implemented to se- 
lect galaxies at more intermediate redshifts (1.5 < z < 3) 
based on optical and infrared (IR) broad-band photometry. 
They include for instance the BzK crit erion (Da ddi et all, 
2004D, the selectio n of BM/BX sources dSteidelet al. , 2004; 

erger et aL| |2004[) and Distant Red Galaxies (DRGs, 
Franx et al.j |20Q3D, the Extremely Red Objects (EROs, 



T hompson et al ^|l999[) or the identification of massive 



sources through the shape of their stellar bump signature in 



the rest-frame near-Infrared ( Simpson fc Eisenhardt , 1999 



Sawicki[ |2002[ [Huang et al.[ |2004p. Each of these meth- 



ods was designed to specifically isolate sub-populations of 
high-redshift galaxies based on a well-defined characteristic 
of their spectral energy distribution (e.g., the blue color of 
their UV continuum, the red color of their evolved stellar 
populations, signatures of star formation reddened by dust, 
...). In this context they have been quite successfully used 
over the past 10 years and their high efficiency has dramat- 
ically revolutionized our understanding of galaxy formation 
in the early Universe. 

The main advantage of these methods is their straight 
forward application which requires only a small number of 
observing bands. Conversely, their calibration over a re- 
stricted wavelength range also implies that they can hardly 
account for the global diversity of galaxy properties at 
high redshift, making perfidious the use of these techniques 
beyond the scope for which they were originally defined. 
In particular it is known that observations at rest-frame 
UV/optical wavelengths can be strongly affected by ex- 
tinction. These techniques could thus suffer from biases 
due to the effect of dust, especially at the bright end of 
the bolometric luminosity function of star-forming galax- 
ies where the contribution of dust-enshrouded star forma- 
tion becomes dominant. Indeed, observations of the deep 
Universe taken with the Spitzer Space Telescope revealed 
the existence of a high-redshift population of dusty lumi- 
nous sources yet almost invisible at shorter wavelength s and 
thus e s caping the traditional UV /optical selections (Houck 
et al I |2Q05| |Dey et al.| |2008| ). What is the fractioiTor 



star formation density produced at high redshift by these 
highly obscured sources? To what extent do galaxies cur- 
rently missed by the standard color selection techniques 
contribute to the growth of stellar mass? 



Over the last decades, observations in the Mid- and 
Far-Infrared as well as at sub-millimeter wavelengths en- 
abled the identification of a large number of luminous and 



large cosmological distances (e.g.. 


Small et al.| 1997 Hughes 


et al.||1998| 


Aussel et al.||1999t|C 


tiary & Elbazl 2001t 


Blain 


et all 12002 


|Le Floc'h et al.| |2004| |Marleau et al.| 


20041 


Coppin et al.| 2008| Capak et al.| 2011). Detailed studies 



spectroscopic surveys have revealed that the bulk of their 
infrared luminosity originates from intense episodes of mas- 
sive star formation, while they can also host powerful active 
nuclei trigg ered by nucelar accretion onto their central black 



holes (e.g., Yan et al.| 


2005| 


Pope et all 2008 


jMenendez- 


IDelmestre et al.| 2 


009 


Hain 


line et al.| |2009| | 


besai et al.| 


2009; Fadda et al.( 


201( 


) ) . Although they have 


Decome rare 



numerous in the past history of the Universe, dominat- 
ing the comoving infrared energy density beyond z ^ 0.5 
and making up to 70% of star-forming activity at z ~ 1 



(Le Floc'h et al. , 2005). At higher redshift their contribu- 



tion could be even larger ( [Caputi et al.[ |2007[ [Rodighiero 



et al. , 2010), although the lack of sensitivity of mid-IR 



and far-IR experiments has prevented reaching definitive 
conclusions on this issue. In fact the respective contribu- 
tions of dust-obscured and unobscured star formation to 
the growth of stru ctures has been a long deb ate for more 
than 10 years (e.g. Adelberger & Steidel, 2000| . Whereas a 
general consensus recognizes that dusty and luminous star- 
forming galaxies played a critical role in driving massive 
galaxy evolution, it is also clear that at the peak of galaxy 
formation the faint-end slope of the UV luminosity function 
was much steeper than observed in the local Universe, im- 
plying a larger contribution of faint galaxi es to the cosmic 
star formation density (Reddy et al., 2008). Unfortunately, 
far-IR observations have not enabled the direct probe of 
these faint high-redshift galaxies yet, while the character- 
ization of luminous star-forming galaxies solely based on 
UV observations requires large extrapolations due to dust 
extinction. So far it has been therefore difficult to reach 
consistent pictures from deep surveys performed at UV and 
far-IR wavelengths. 

The goal of this paper is to address the fraction of lu- 
minous high-redshift galaxies that may be missed by the 
different color selection techniques commonly used to iden- 
tify distant sources, and to quantify the impact of this bias 
by estimating the contribution of these objects to the IR lu- 
minosity and cosmic star formation density at 1.5 < z < 3. 
We carried out this work using the deep 24/im qbservar 



tions (F24um > 0.08 mJy) of the COSMOS field dScoville 
et al.| | 2007b[ ) obtained with the MIPS instrument ( [Kieke 



et al.| |2004p on-board the Spitzer Space Telescope. Among 
all facilities allowing the probe of dusty galaxies in the 
distant Universe (i.e. SCUBA, AzTEC, LABOCA, MIPS 
70 and 160/im, ...), MIPS 24/im observations provide one 
of the deepest sensitivity limits currently reachable up to 
z ~ 2 — 3. Furthermore the deep optical/near-IR imaging as 
well as the high-qual ity photometric r edshifts obtained in 



the COSMOS field ( ,Ilbert et al. 



characterization of the redshift distribution associated to 



2009|) enabled a detailed 



These sources are referred as Luminous Infrared Galaxies 
(LIRGs: lO^^Lo < LiR < lO^^L© with Lir = Lg -loooum) and 
Ultra-Luminous Infrared Galaxies (ULIRGs: Lir > 10 ^L©). 
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the whole s ample of 24/im sources up to z ~ 3 (Le Floc'h 
et al. , 2009), hence providing a complete and unbiased view 
of the population of dusty high-redshift sources selected 
at F24/im > 0.08 mJy. Taking advantage of the large cov- 
erage of COSMOS to minimize the effect of cosmic vari- 
ance, we applied the color criteria associated to four differ- 
ent selections of distant galaxies {BzK & BM/BX sources, 
''IRAC peakers'' and ''Optically- Faint IR-hrighV' (OFIR) 
objects), and we quantified the amount of high-redshift 
24/im sources missed by each of these techniques. Our data 
are described in Sect.[2j the criteria that we explored are 
presented in Sect. [3] and the corresponding sub- selections 
that we obtained based on the MlPS-detected COSMOS 
population are discussed in Sect. 4. In Sect. |5] we present the 
rest-frame 8/im luminosity function (LF) of the complete 
high-redshift MIPS sample as well as the contribution of 
the sub-populations of galaxies respectively identified with 
the color selection techniques mentioned earlier. Our results 
are discussed in Sect. [6] and we finally present our conclu- 
sions in Sect.[7[ Throughout this paper we assume a ACDM 
cosmology with Hq = 70 km s~^, Vtrn = 0.3 and = 0.7. 
Magnitudes are given in the AB system and the conversions 
between luminosities and Star Formation Rates are com- 



29 395 sources detected at 24/im over an effective surface of 
1.68 deg^ 

The counterparts of the MlPS-selected detections were 
identifie d following the sam e procedure as the one described 
by |Le Floc'h et al. (2009). Given the much higher den- 
s ity of sources dete cted at optical wavelengths in COSMOS 
(Capak et al. , 2007) compared to those detected with MIPS, 
a direct cross-correlation between the 24/im-selected cata- 
log and the optical observations could lead to a large num- 
ber of spurious associations with optically-detected galaxies 
randomly aligned close to the line of sight of the MIPS 
sources. Hence, we first cross-correlated the 24/im data 



puted assuming the Initial Mass Function from |Salpeter 
(IT9551). 



2. The Data 

The sample of luminous star-forming galaxies used in this 
work was selected from the deep S pitzer /MIPS observa- 



tions of the 2 deg^ COSMOS field (jSanders et al. , 2007) 



and it is based o n the 24/im source catalog obtained by 



Le Floc'h et al.| ([2009|. COSMOS is the largest contigu- 
ous im aging survey ever undertaken with the HST ( Scoville 



et al. 



2007a} to a depth of i+ =28.0 [AB] (5 ( 7 on a n op- 



timally extracted point source; Capak et al. 2007^. The 

the 



COSMOS field is located near the celestial equator to en- 
sure visibility by all ground and spaced-based astronom- 
ical facilities. It is devoid of bright X-ray, UV, and ra- 
dio foreground sources, and compared to other equato- 
rial fields the Galactic extinction is remarkably low and 
uniform {{E(^b-v)) — 0.02). Extensive multi- wavelength 
follow-ups with ground-based and spaced-based facilities 
have also been achieved with high-sensitivity imaging and 
spectroscopy spa nning the entire spectrum from the X- 



ray to the radio ( Hasinger et al.[ 2007; Schinnerer et al. 



2007 
deep 

ent b ands 

(e.g., [Taniguchi et al.[ |2007[ [Capak et al.[ |2007[) 



[Lilly et al.[ |2007| [Elvis et al.[ |2009). In particular 
UV-to-8/im photometry was obtained over 30 differ- 



using 



narrow, medium and broad-band 



filters 
while 

far-IR and submillimeter /millimeter observations were per- 
formed w ith various inst r ument s such as MIPS on-board 
Spitzer ([Sanders et al.| [ 2007| ), PACS and SPIRE on- 

al. submitted, Oliver et al. in 



(Sanders et al. 
board Herschel (Lutz et 



prep.), SCUBA-2, AzTEC, MAMBO and BOL OCAM (e.g. 



Bertoldi et al., 2007 Austermann et al. 2009). 



While the 24/im source catalog of COSMOS covers a 
total area of 2deg^, we restricted our sample to the field 
outside the regions contaminated by very bright or satu- 
rated objects and where the photometry is less accurate. We 
also applied a conservative flux cut of 80/iJy, corresponding 
to a compl eteness of ^ 90% in t he source extraction per- 
formed by |Le Floc'h et al.| ( [2009| ). This led to a sample of 



with t he COSMOS i^^-band catalogue of McCracken et al. 
(12010 ) to minimize the risk of wrong associations. The den- 
sity of sources in this catalogue is substantially smaller 
than at optical wavelengths, while it is deep enough (5<j 
for i^s=23.7) to allow the identification of near-IR counter- 
parts for most of the 24/im detections in COSMOS (e.g., 
Le Floc'h et al.[[2009 ). Also, the i^^-band observations were 
carried out under very good seeing conditions (^0.7^' at 
2.2/im), leading to a PSF much narrower than obtained in 
the other near-IR images of COSMOS (e.g., IRAC-3.6/im) 
and allowing more robust identifications in the case of 
blended sources. Our correlation between the 24/im and 
Kg-band data was performed with a matching radius of 2". 
Given the width of the 24/im PSF (FWMH-6"), this radius 
allowed us to identify most of the 24/im sources while also 
minimizing the nu mber of multiple matches. Similar to the 
results obtained by Le Floc'h et al. (2009), only 765 sources 
from the initial sample of 29 395 mid-IR objects could not 
be matched in the near-IR, while up to 84% sources from 
the MIPS catalogue were identified with a clean and single 
Kg-band counterpart. For the rest of the sample (13%), two 
possible matches were found within the matching distance. 
For those cases we decided to keep the closest counterpart, 
as for most of them (~60%) the centroid was at least twice 
closer to the 24/im source than the second possible match. 

In a second step we correlated the list of these Ks- 
band identifications with a recently-updated version of the 
-band selected c atalogue of photometric redshifts from 
Ilbert et al.[ (|2009[) using a matching radius of V . Given 
the depth of the C 



OSMOS i+-band data (e.g., [Capak et aO 
2007^ at the end 1306 MIPS sources could not be matched 
at optical wavelengths, leading to the identification of opti- 
cal counterparts and redshift determination for more than 
95% of the initial 24/im-selected sources. As part of this 
second cross-correlation, double optical matches were ob- 
tained for only ~1% of the K^-band counterparts, under- 
lying the robustness of our identification at these short 
wavelengths. For these few cases we kept ag ain the clos- 
est possible optical association. The catalog of jllbert et al. 
(2009) includes the band- merged COSMO S photometry 
published by jCapak et al.| ( |2007D , |McCracken et al.[ ( [2010 ) 
and Jlbert et al.j ( |2010p from the (/-band to the IRXC 
8/im. Our cross-correlation thus yielded direct determina- 
tion of optical/near-IR broad-band magnitudes for all se- 
lected galaxies. In particular, we considered the photome- 
try in the li*, 5j, Vj, r+, i~^ and z~^ band^ where the 5-<j 
magnitude limi t reaches respectively 26.5, 26. 6, 26.5, 26.6 



26.1 and 25.1 ( [Capak et aT[ [2007] ). Finally, [Ilbert et al" 



In this work we also applied to the COSMOS broad-band 



photometry the systematic offsets inferred by [Ilbert et al. ( |2009 ) 
when computing their photometric redshifts (see their Table 1). 
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( 2QQ9| ) also provide the identific ation of sources detected at the CFHT and at the Subaru Telescope with the u*^ 
witF XMM oyer CO SMOS (e.g., |Brusa et al.] |2Q1Q[ |2QQ7| and r+ filters (Capak et al] |2QQ7|), which d iffer substan- 



Salvato et al. 2009), which allowed us to systematically ex- tially from the filters used by Steidel et al. (12004). As a 



elude the X-ray detected AGNs from our final sample down 
to a flux limited S'o.5-2feey 
than 4%). 



5 X 10 ^ erg cm ^ s (less 



3. Photometric selection techniques 

3.1. Optical and NIR color selections 
3.1.1. The BzK selection 

Based on deep phot ometry obtained wit h B-band, 2;-band 
and i^-band filters, [Daddi et al. ( |2004 ) proposed a crite- 
rion to select star-forming galaxies in the redshift range of 
lA<z < 2.5: BzK = {z-K)ab-{B-z)ab > -0.2. In the 
plane defined by 2; — as a function of B — z (hereafter the 
BzK diagram), the extinction vector obtained for the at- 
tenuation law of Calzetti et al. (2000| is parallel to the line 
characterized by BzK = —0.2 at 1.5 < z < 2.5. Therefore 
this criterion presents the strength of being mostly inde- 
pendent of dust obscuration. 

The 5-band^-band and i^^-band photometry used by 



Daddi et al. 



( |2004| ) was based on the Bessel-5 and Kg fil- 

ters available at the VLT as well as on the F850LP z-band 
filter of HST. On the other hand, our COSMOS photome- 
try was determined through Bj and 2;+-band observations 
obtained a t the Subaru Telescope and Kg-band data ta ken 
at CFHT (ICapak et al] |2007| [McCracken et al.| |2010l). In 
order to adapt the original BzK criterion to the COSMOS 
data, we thus determined the evolution of the BzK vari- 
able through the COSMOS filters as a function of redshift 
and with a set of star-formin g galaxy templates similar to 



the library of SEDs used by Ilbert et al. (2009) for their 



catalog of photometric redshifts in COSMOS, and we com- 
pa red this BzK variable to the BzKq color initially defined 
by iDaddi et al.| ^20041 ). The star-forming SEDs used by 
Ilbert et al. ( |2009 ) wer e generated with th e stellar popula- 
tion synthesis model of'Bruzual fc Char lot | ([2003) assuming 
exponentially-declining star formation histories and star- 
burst ages between 0.03 and 3 Gyr. In the redshift range of 
lA < z < 2.5 we found BzK — BzKq = 0.1 mag with a very 
small dispersion driven by the choice of SED template. To 
apply the BzK selection technique to the COSMOS sur- 
vey we thus added a system atic offset of +0.1 m ag to the 
initial criterion proposed by Daddi et al. ( [2004 ), leading 
to the identification of star-forming BzK galaxies with the 
following color: 



BzK > -0.1 



3.1.2. The BM/BX criteria 

Another widely-used select ion of distant 
sources was proposed by 'Steidel et al 



(1) 



s tar- fo rming 
(20041) and 



[Adelberger et al.|(|2004) using UV/optical broad- band color 



selection techniques. Based on photometry obtained in the 
Un^ G and R bands and using the associated Un — G and 
G — R colors, they defined two criteria refered as the "BM" 
and "BX" selections to identify galaxies at 1.4 < z < 2.1 
and 1.9 < z < 2.7 respectively. 

At the wavelengths where these two selection techniques 
can be applied, the COSMOS observations were carried out 



result we were not able to adapt the initial BM/BX selec- 
tion to the COSMOS photometry by applying simple terms 
of color corrections and/or systemat ic offsets to ou r data . 
Following the approach proposed by Grazian et al. (2007) 
we considered a slightly modified version of the BM/BX 
criterion using the ^i* — Vj and Vj — i+ co l ors ins tead of the 
colors initially used by Adelberger et al. (2004). To deter- 
mine the most appropriate color cuts for our it*Vjz+' dia- 
gram, we computed the evolution of color tracks followed 
by star-forming galaxies as a func tion of redshift, in a way 
very similar to that employed by [Adelberger et al. ( 2004 ) 
so as to ensure a selection of star-forming sources as close 
as possible to the original criterion. We used the spectral 
energy distributions of two starburst galaxies among the 
list of galaxy templates considered by Ilbert et al. (2009) 
for their analysis of the COSMOS photometric redshifts. 
These two SEDs were chosen as being the most repre- 
sentative templates of galaxies beyond z ^ 1 according 
to the SED fitting infered by |Ilbert et all (I2009|). These 
two templates correspond to starburst models taken from 
dBruzual fc Chariot ||2QQ3l hereafter BC03) to which emis- 
sion lines has been added to have a better understanding 
of th e COSMOS colors (for more details see Ilbert et al. 



2009). Following Adelberger et al. (2004), we also took into 
account the effect of extinction by reddening each template 
up to E{B — V) = .2 assuming the attenuation curve of 



[Calz etti et all ( |2000D . 

The exact determination of our new selection is illus- 
trated on Fig[l] First we defined the lower and upper limits 
of our criterion on the ^x* — Vj color using the location of 
the galaxy tracks at respectively z = l.b and z = 2.8: 



ii* - Vj > 0.52(V> -0.1 
ii*-Vj<2.5(Vj-i+) + 1.0 

Second, we used a lower limit of —0.25 for the Vj 



(2) 



color. It is similar to the cut chosen by Grazian et al. (2007), 
although we note that its effect is negligible since very few 
sources are bluer than this limit. Finally, we imposed an 
upper limit on this second color so as to allow the selection 
of sources corresponding to the reddest templates of our 
diagram {E{B — V)=0.2mag): 

ii* - Vj < -12.5(Vj - z+) -h 9.0 (3) 

Sim ilar to what had been noticed by [Adelberger et aT] 
(12004), the galaxy tracks illustrated in Fig [l] reveal an 
important degeneracy between the colors of low-redshift 
galaxies and more distant sources affected by dust extinc- 
tion. The Vj — upper limit that we adopted corresponds 
to the compromise of allowing the selection of dusty galax- 
ies while minimizing the contamination of sources at z < 1. 
As we will see later in this Section, this will have strong 
impact on the selection of luminous and dusty star-forming 
objects with E{B-V) > 0.2 using the BM/BX technique. 
We kept however a strict cut at E{B — V) = 0.2 not to al- 



ter th e original BM/BX criterion proposed by Adelberger 



et al. 



( [2004D 



To verify the robustness of our new color selection we 
overploted in Fig [l] the observed colors of the z+-band se- 
lected COSMOS galaxies with the aforementioned redshift 
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Fig. 1. Selection of galaxies at 1.5 < 2; < 2.8 based on li* — Vj and Vj — z+ colors. The galaxy tracks represented by 
the green and red solid curves correspond to the colors predicted as a function of redshift for two typical star-forming 
galaxy templates used by Ilbert et al.| ( 2009) for their determination of photometric redshifts in COSMOS. The dashed 
an d dotted lines s how t he effect of dust extinction produced by respectively E{B — V) = 0.1 and 0.2 assuming the law 
of Calzetti et al. (2000). The cross symbols along the tracks correspond to z = 0, 1, 1.5, 2, 2.8 and 3. The black solid 
lines show the limits of our modified BM/BX criterion. The distribution of the optically-selected COSMOS galaxies at 
1.5 < z < 2.8 (black dots, limited for clarity to sources with i+ < 25 and photometric uncertainties smaller than 0.1 mag 
in li*, Vj and i+) illustrates the high reliability of our color selection. 



cut of z> 1.5 dCapak et aL| |2007[ [Ilbert et al.l|2009D . For 

clarity, only sources with < 25 are shown, so as to mini- 
mize the dispersion coming from photometric uncertainties. 
As expected from the SED tracks, we see that our crite- 
rion is globally well suited to the identification of optically- 
selected star- forming galaxies at 1.5 < z < 2.8. Out of the 
28 141 i+-band sources located in this redshift range and 
selected at i+ < 25 from the photometric redshift catalog 
of Ilbert et aL ( 2009) , 24 826 sources (i.e., 88%) satisfy our 
modified BM/BX selection. Among the remaining objects 
that escape the selection, some are substantially affected by 
dust extinction and their optical colors are thus redder than 
allowed by the criterion, while some others are character- 
ized by photometric uncertainties moving them just outside 
of our BM/BX selection box. Most of them, however, corre- 
spond to blue sources with photometric redshifts just above 
z = 1.5. Therefore the small incompleteness of our criterion 
mostly reflects the color degeneracy that appears at this 
redshift because of SED variations and dust extinction, as 
well as uncertainties affecting our photometric redshifts in 
the COSMOS field. Relaxing our lower limit on the ii* — Vj 



to reduce this incompleteness of 12% would inevitably in- 
crease the fraction of contaminants from sources at z < 1.5. 



3.1.3. Distant "stellar bump dominated" galaxies 

After the successful launch of the Spitzer Space Telescope, 
the advent of deep surveys carried out with the IRAC in- 
strument enabled a new approach for identifying distant 
star-forming galaxies, based on the rest-frame 1.6/im bump 
produced by the H~ opacity minimum in the atmosphere s 



Sawicki, 2002). 



Dump is shifted 



of cool stars (Simpson & Eisenhardt, 1999 
At 1.5 < z < 3 the signature of this stellar 
between 3 and 9/im and it reveals itself through specific 
colors across the IRAC bands. It allows the identification 
of distant galaxies in a way that is less subject to extinc- 
tion compared to the optically-based selections discussed 
earlier. 

Several selection techniques have been explored to apply 
this method to large galaxy samples. They involve either a 
direct fit of the IRAC fluxes with stellar templates (e.g., 
Sorba & Sawicki, 201 0]), the use of two-colo r criteria (e.g., 
iPapovich^ i2008, ,Huang et al.[ |2004[ |2009D or the simple 
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identi fication of the SEP pe ak across the 4 IRAC channels 



(e.g., Lonsdale et al., 2009). For our current analysis we 



considered the selection introduced by Huang et al. (2004). 
It is based on the color cuts 0.05 < [3.6J - [4.5J < 0.4 
and -0.7 < [3.6] - [8.0] < 0.5, where [3.6] denotes the 
AB magnitude in the 3.6/im band (and likewise for the 4.5 
and 8.0 /im bands). Extensive follow-up have confirmed the 
re liability of this criterion for selecting g alaxies at 1.5 < z < 
3 dHuang et al] |2009l [Pesai et al.j|20Q9| . These sources will 
b e refer ed as the IRAC Peakers hereaft er . 



3.2. Selection of Optically Faint IR-bright sources at z^2 

While the selection techniques discussed in the previous sec- 
tion pertain by definition to sources primarily identified at 
optical and near-IR wavelengths, observations undertaken 
in the thermal Infrared and the submillimeter have revealed 
a large number of luminous galaxies associated with optical 
counterparts much fainter than the typical magnitude lim- 
its considered in the surveys of the distant U niverse c arried 
out in the visible (/ ~ 25 mag AB, e .g., Hughes e t al.[ il998; 



Houck et al.[|2005[ryan et al.[|2005| ). Given the extremeTH"- 
to-optical colors that result from these faint magnitudes 
at visible wavelengths, various criteria based on the mid- 
IR/optical flux ratio have thus been proposed for selecting 
dusty star-forming sources at 2; > 1. However the relevance 
of the galaxy sub-samples selected with this approach and 
their contribution with respect to the global population of 
high redshift sources has barely been explored so far. Here 
we describe the selection technique that was proposed by 
Dey et al. ( 2008 ) , based on the ratio between the 24/im and 
the optical i^-band flux densities. Using the Spitzer obser- 
vations of the Bootes Field and the NOAO Deep Wide- 
Field Survey (NDWFS) they analyzed the population of 
24/im sources brighter than 0.3 mJy (i.e., the depth of the 
MIPS imaging in Bootes) and they deflned a selection of 
dust-enshrouded high redshift source^ with the following 
criterion: 



R 



(4) 



[24] ^ 7.4 mag (AB) 

where [24] refers to the AB magnitude measured in 
the MIPS-24/im bancQ Using spectroscopic follow-up per- 
formed with the IRS spectrograph on-board Spitzer and 
with the LRIS/DEIMOS instruments at the Keck tele- 
scopes, they have shown that their technique provides a 
reliable selection of dusty luminous galaxies at 1.5 < z < 3, 
with a small fraction of contaminants from lower redshifts. 

4. Application to the 24/xm selected galaxy 
populations in COSMOS 

4.1. The optical/near-IR selections applied to the 24jj.m 
galaxy population: associated redshift distributions 

To quantify how the UV/optical selections of distant 
sources suffer from biases and incompleteness due to dust 
extinction, we applied the various criteria previously dis- 
cussed to the MIPS-24/im sample described in Sect.|2] We 

^ refered as "Dust-Obscured Galaxies" in their selection at 
F24fxm > 0.3 mJy. 

The relation given by "Dey et al. (^20081 uses the Vega sys- 
tem and it was converted following the MIPS Data Handbook: 

miVega) = -2.5 x /o^io^^^f#^. 



then analyzed the relative weight of the different sub- 
samples of 24/im sources identified with the BzK^ BM/BX 
and IRAC peaker selections with respect to the whole pop- 
ulation of 24/im-selected galaxies. 

Our results are illustrated in the different panels of 
Fig. 2. On one hand, Fig.|2^ shows the redshift distribu- 
tions associated with our three UV/optical selections and 
compared to the redshift distribution obtained for the en- 
tire 24/im source population in COSMOS. On the other 
hand, Figs.[2]3, |2]3 & |2]i represent the corresponding color 
diagrams along with the distribution of the 24/im sources 
independently identified in the redshift range associated 
to each selection (1.4 < z < 2.5 for the BzK selection, 

1.4 < z < 2.8 for our modified BM/BX criterion, and 

1.5 < z < 3 for the IRAC pea kers) using our reds hift iden- 
tifications from the catalog oflll bert et al.| (|2QQ9[). The un- 
certainties characterizing the photometric redshifts of the 
24/im galax y population in COSMO S h ave been thoroughl y 



discussed by |Le Floc^h et al.| ( |2QQ9| ) and |Ilbert et aL| ( |2QQ9 ). 

Even at 1 < z < 3 they are substantially smaller than the 
typical redshift range where the color selections operate 
{l(7z < 0.15 up to < 25 mag and z < 3), and therefore 
they should not affect our global conclusions. Not surpris- 
ingly we see that the distributions characterizing the BzK, 
BM/BX and IRAC peaker sub-samples globally cover the 
redshift range of 1 < z < 3, as expected from the definition 
of their associated selections. Nonetheless, these distribu- 
tions clearly show that some criteria present substantial 
biases. Hereafter we discuss in more details the redshift 
properties characterizing each individual sub-sample. 

• BzK galaxies: Out of the 7227 MIPS sources lying in 
the redshift range 1.4 < z < 2.5, we found that 6 623 
objects satisfy the BzK criterion. This represents a 
fraction of 92% of the whole MIPS population at these 
redshifts, implying that the BzK selection is particu- 
larly efficient and weakly affected by the effect of dust 
obscuration. This high efficiency can be explained by 
the fact that in the BzK diagram the extinction vector 
evolves parallel to the line defining the BzK criterion 
for star- forming galaxies {Bzk = —0.1), hence allowing 
the selection of deeply obscured sources provided 
they are detected in the i^T-band. Furthermore, the 
redshift distribution (Fig. 2a) reveals that the presence 
of contaminants at z < 1.4 is clearly negligible, while 
the BzK colors of the MIPS sources identified in the 
lA < z < 2.5 redshift range (based on photometric red- 
shifts) but not selected as BzK sources are distributed 
very close below the threshold of BzK = —0.1 (see 
Fig. 2b). It shows that the small incompleteness of the 
BzK selection at 1.4 < z < 2.5 could simply result from 
optical/near-IR photometric uncertainties spreading 
sources below the borderline of the criterion, as well 
as small uncertainties on the photometric redshifts of 
sources close to the redshift boundaries of the selection. 
As an interesting aside we also note that 84% of the 
MIPS sources with 2.5 < z < 3 are still selected as 
BzK sources. It suggests that the BzK criterion could 
also be applied successfully at even higher redshift, at 
least for our sample of dusty luminous galaxies. 
It is important to stress that the high success rate of 
identifications that we find with the BzK technique 
is largely due to the depth of the optical and near-IR 
observations available in the COSMOS field. A large 
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Fig. 2. (a) Redshift distributions of the MIPS galaxies identified with the BzK (red), BM/BX (green) and IRAC 
peaker (orange) selection techniques, compared to the distribution obtained for the whole sample of COSMOS sources 
with F24^rn > 0.08 mjy (black solid histogram), (b) The BzK color distribution of the MIPS galaxies identified at 
lA < z < 2.5 with the COSMOS photometric redshifts. The BzK star-forming sources (S-BzK) are identified above the 
diagonal solid line, (c) Distribution of the 1.5 < z < 2.8 MIPS galaxies in the Vj — versus — Vj color diagram. 
The area defined by the solid line corresponds to our modified BM/BX criterion (see Sect. 3.1.2). (d) IRAC colors of the 
MIPS sources located at 1.5 < z < 3. The box illustrates the IRAC peaker selection criterion proposed by [Huang et al. 
(2004). The typical uncertainties for each population are shown in the corners of the figures. 



number of high redshift 24/im galaxies are indeed 
associated with faint optical/near-IR counterparts and 
the ratio of MIPS sources identified with the BzK 
selection would have been much lower if we had used 
shallower data at optical and near-IR wavelengths. 
As an example we restricted the COSMOS i^g-band 
sample to respectively Ks =22, Ks =22.5 and Ks =23, 
and we found that with these additional cuts only 42%, 
72% and 89% of the MIPS population at 1.4 < ^ < 2.5 
would have been identified with the BzK selection. 
Given the large contribution of dusty luminous galaxies 
to the cosmic star formation density at z ~ 2 (e.g., 
Caputi et aT] |2QQ7[ [Rodighiero et al.| |2Q1QD this clearly 
illustrates the critical need for deep near-IR obser- 
vations when probing the star-forming high-redshift 
galaxy population with the BzK criterion. 
Finally, we also recall that 4% of sources from the very 
first sample of 24/im detections extracted from the 
MIPS imaging of COSMOS (see Sect. |2| could not be 



matched with the optical catalog of Ilbert et al. (2009). 
We then have no information on their colors and their 
redshift. Given their faint ness at short wavelengths 
these galaxies most probably lie at z>l. If they were all 
located at 1.4 < z < 2.5 they would represent 15% of 
the whole population of 24/im sources identified in this 
redshift range. The fraction of MIPS sources selected 
as BzK galaxies would thus decrease from 92% down 
to -78%. 

BM/BX sources: The distribution shown in Fig [2k 
reveals that the BM/BX criterion enables the identm- 
cation of dusty galaxies in the redshift range expected 
from the analysis that we presented in Sect. 3.2, but 
with a much lower efficiency than the BzK and the 
IRAC peaker selections. Out the 7459 MIPS sources 
with 1.5 < z < 2.8, only 3 754 (i.e., --50%) are indeed 
identified as BM/BX galaxies. To understand the origin 
of this bias we illustrate in Fig. [2b the distribution 
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of the 24/im sources at 1.5 < z < 2.8 in the BM/BX 
diagram. It clearly shows that a large fraction of them 
is located outside of the BM/BX selection area that 
we defined earlier, because of much redder colors than 
observed on average in the optical sample selected at 
1.5 < z < 2.8. Over this redshift range, the comparison 
between the BM/BX properties of the MIPS sources 
and that of the COSMOS i+-band selected population 
displayed in Fig.[T] is in fact particularly striking. 
According to the galaxy tracks represented in this 
figure, the MIPS sources located to the right hand 
side of the BM/BX selection box correspond either to 
galaxies with substantial extinction {E{B — V) > 0.2) 
or to lower redshift contaminants that we could have 
falsely identified with high-redshift objects. The later 
is however unlikely given the very low fraction of 
catastrophic failures among the photometric redshifts 
of the MlPS-selected population i n COSMOS (Ilbert 



et al. , 2009 



2009), as well as the ex- 



Le Floc'h et al. , 

cellent agreement that we just found between the BzK 
selection and the MIPS sample at z ~ 2. Therefore we 
conclude that the large fraction of 24/im sources missed 
by the BM/BX selection correspond to galaxies truly 
located at 1.5 < z < 2.8 but strongly reddenned by 
dust extinction. This is obviously not surprising since 
the Spitzer mid-IR observations naturally favor the 
selection of dusty objects. Yet it clearly reveals how the 
BM/BX criterion is biased toward the identification of 
UV-bright sources with no or little dust obscuration 
{E{B — V)< 0.2) and how this selection can thus miss 
a large amount of luminous dust-enshrouded galaxies in 
the distant Universe. In Sect. 5 we will quantify in more 
details the impact of this bias on the bright end of the 
IR luminosity function of star- forming galaxies at z ~ 2. 



• IRAC peakers: Among the 7 755 MIPS sources located 
at 1.5 < z < 3 a total of 4 942 galaxies were se- 
lected as IRAC peakers. This corresponds to an average 
completeness of ^64% over this redshift range, show- 
ing that this criterion is also suited to distant galaxies 
but slightly less efficient than the BzK selection. This 
smaller efficiency must be due to other effects than dust 
extinction since the IRAC criterion is based on data 
taken at longer wavelengths. A possible explanation is 
that the color selection of IRAC peakers has been re- 
stricted to identify only the galaxies where the stellar 
bump is particularly pronounced, while dusty objects 
usually show a wider diversity of SEDs in the near-IR 
given the combined effect of stellar and hot dust emis- 
sion. Besides, the COSMOS photometric uncertainties 
in the IRAC bands are typically la rger than those ob - 



tained at shorter wavelengths (e.g., Ilbert et al. , [2009 ) 
As we can see in Fig. 2d this may result in distributing 
galaxies at 1.5 < z < 3 over a slightly wider range of 
IRAC colors while increasing the contamination from 
galaxies at lower redshifts. Indeed we note that the red- 
shift distribution of the 24/im sources identified as IRAC 
peakers extends below z ^ 1.5 with 23% of the sub- 
sample being located at 1 < z < 1.5. Finally, although 
bright X-ray AGNs were removed from our initial sam- 
ple, we can not fully exclude a residual contamination 
from high-redshift obscured nuclei detected at 24/im. 
Their SED is usually characterized by a rising feature- 
less hot dust continuum, which can totally outshine the 



stellar e mission of the galaxy when the AGN is luminous 
enoug h (|Bran d et al." 2006 Menendez-Delmestre e t ah] 
2009j Desai et al., 2009). In this case their host exhibit 
much redder colors than those allowed by the "IRAC 
peaker" selection, which would contribute to lowering its 
efficiency. However, the AGN contribution to the whole 
population of IR galaxies and to the Cosmic I nfrared 
Background is likely not larger than 15% (e.g., Jauzac 



et al. , 2011 ). Hence, their possible contamination in our 



sample can not be the main explanation for the incom- 
pleteness of the IRAC peaker method. 

4.2. Infrared-luminous sources with optically-faint 
counterparts at z^2 

In the analysis presented in the previous sub-section we 
have seen that a substantial fraction of the high-redshift 
mid-IR selected galaxies are strongly reddened by dust ex- 
tinction, and their identification can thus be missed when 
purely relying on rest-frame UV color selection techniques. 
Conversely, the selection of distant sources characterized 
by a large excess of IR emission such as the Optically- Faint 
IR-bright (OFIR) objects will necessarily be biased against 
dust-free galaxies. To quantify this effect, we applied the 
optical/mid-IR criterion described in Sect. 3.2 to our MIPS 
galaxy sample using the r+-band observations of COSMOS 
and we extended this criterion down to F24^m = 0.08 mJy 
in order to probe a larger range of mid-IR luminosities. In 
Fig. pi we show the distribution of the r+ — [24] colors for 
the MlPS-selected sources as a function of redshift, which 
confirms that most of the Op tically-Fain t IR-b right galax- 
ies satisfying the criterion of Dey et al. (2008) are indeed 
located at 1.5 < z < 3. Figs.[3]3 illustrate the corre- 
sponding redshift distributions obtained for MIPS sources 
with respectively F24,^rn > 0.08 mJy and F24/im > 0.3 mJy, 
compared to the distribution of the full MlPS-selected sam- 
ple in COSMOS. The Optically-Faint IR-bright sources 
only represent ~30% (2 350 sources on the 7 755 sources 
with 1.5 < z < 3) of the population of MIPS galaxies se- 
lected with F2Anm > 0.08 mJy at 1.5 < z < 3, but as we 
will see in the next section their contribution is rising with 
the mid-IR luminosity. Above a 24/im flux of 0.3 mJy their 
fraction reaches ~60% of the 24/im-selected population at 
^ - 2. 

The fact that the fraction of galaxies with 24/im to 
optical-band flux ratios satisfying Eq. [4] rises as a function 
of 24/im flux is not an artiflcial selection effect due to the 
relative depths of our MIPS and r+-band COSMOS obser- 
vations. Even at the 24/im flux limit of 0.08 mJy the optical 
COSMOS data are deep enough to identify galaxies as red 
as the color thr eshold used in the criterion proposed by 
Dey et al.| (|2008|. Also it implies that the 4% of the initial 
24/im detections that we could not identify at optical and 
near-IR wavelengths (Sect. 2) satisfy Eq.|4] independently of 
their mid-IR flux density. Since these extremely red colors 
are unlikely associated to galaxies at z < 1 (see Fig.[3|i) 
the relative contribution of galaxies selected at 1.5 < z < 3 
with this criterion could reach ~45% at F24^m > 0.08 mJy. 

As we already argued this selection of high-redshift 
galaxies based on extremely red mid-IR/optical colors is 
by construction biased toward dust-obscured sources with 
very faint optical luminosities. For example we found that 
66% of these objects have r+ > 26 mag while the ri- 
band magnitude distribution of the MIPS sources with 
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Fig. 3. (a) The r+ - [24] color distribution of the MIPS 
24/im sources as a function of redshift in the COSMOS field. 
Galaxies selected as Optically- Faint IR-bright sources are 
characterized by r+ — [24] colors redder than the thresh- 
old indicated by the dashed line. This criterion effectively 
selects galaxies at 1.5 < z < 3. The typical uncertainties 
on the r+ — [24] color have been calculated at 1.5 < z < 3 
and are shown on the right corner, (b) Redshift distribu- 
tion of the Optically- Faint IR-bright galaxies with F24^rn > 
0.08 mJy (blue filled histogram) compared to the redshift 
distribution obtained for the full MlPS-selected galaxy pop- 
ulation in COSMOS (black solid line), (c) Redshift dis- 
tributions derived from the MIPS sample restricted to 
sources with F2A^m > 0.3 mJy (same colors definitions than 
(b)). The relative contribution of the OFIR galaxies at 
1.5 < z < 3 increases with mid-IR luminosities. 



1.5 < z < 2.8 peaks at r+ ~ 25 mag (only 21% of 
them have r+ > 26 mag). Although they do not repre- 
sent a dominant population in terms of source density, 
this nicely illustrates the contribution of luminous high- 
redshift star-forming galaxies that may be missed by op- 
tical surveys because of dust obscuration. In fact, 65% of 
the Optically-Faint IR-bright sources have an extinction 
of E(B - V) ^ 0.4 a ccording to the optical SED fitting 



of |Ilbert et ah] ( 2009[ ), while this percentage decreases to 
only 44%) for the whole population of MIPS sources in the 
same redshift range. Furthermore, galaxies selected with 
this technique differ dramatically from those identified with 
the BM/BX selection, which we found to be strongly biased 
against dusty high-redshift galaxies. In Fig. [4] we compare 
in the BzK diagram the distribution of the Optically- Faint 
IR-bright sources with the distribution of galaxies selected 
with our modified BM/BX criterion. The overlap between 
the two sub-samples is relatively small as only 8% of the 
MIPS sources at 1.5 < z < 2.8 satisfy the two selection 
criteria, and not surprisingly the Optically-Faint IR-bright 
objects are characterized by much redder colors than the 
BM/BX sources. Since the Optically Faint IR-bright galax- 
ies and the BM/BX selected sources share very similar red- 
shift distributions (compare Figs. 2a & 3b) and given the 
extinction vector at z ~ 2 in the BzK diagram (Fig.[4|, 
the lack of overlap between these two populations is mostly 
due to the effect of dust extinction. 

In Fig.|5] we illustrate in a more quantitative way how 
the Optically Faint IR-bright objects and the other selec- 
tions overlap with each other. Since we found that the high 
efficiency of the BzK criterion is maintained up to z ~ 3 (at 
least when applied to our mid-IR sample, see Sect. 4.1), we 
compared our different selections over the largest possible 
redshift range (i.e., 1.5 < z < 2.8) to minimize the sta- 
tistical uncertainties. More than half of the population of 
Optically Faint IR-bright sources are also selected as IRAC 
Peakers indicating that an important part of the sources 
presenting extreme MIR color excess have an SED domi- 
nated by the stellar bump. 



5. Galaxy S/xm luminosity function and Star 
Formation Rate Density at z ^ 2 

In the two previous sections we described a number of color 
selection techniques that have been widely used for identify- 
ing star- forming galaxies at 1.5 < z < 3, and we quantified 
their respective contribution to the total number density 
of MlPS-selected high-redshift sources in order to estimate 
the bias that these selections suffer because of dust extinc- 
tion. We want to extend this analysis by measuring the 
mid-IR luminosity function of the different sub-samples of 
galaxies selected based on these techniques, so as to infer 
their contribution to the total IR luminosity density of the 
Universe and the cosmic star formation rate density ob- 
served at z ~ 2. Here we detail the methods used to obtain 
the 8/im rest-frame Luminosity Function (LF), the IR lu- 
minosity density and the Star Formation Rate (SFR). 

5.1. Sjim luminosity function 

The galaxy luminosity functions were all computed at rest- 
frame 8/im. The redshift distributions of the sources se- 
lected with the techniques described earlier mostly peak at 
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Fig. 5. Overlap between the populations of 24 jam sources with F24^m > 0.08 mJy on the redshift range 1.5 < z < 2.8. 
The percentages in bold are related to the total number of sources selected by the criterion considered, the percentages in 
italic correspond to the overlap between the populations. Left: Overlap between the BzK, the BM/BX and the Optically 
Faint IR-bright sources. Right: Overlap between the BzK, the Optically Faint IR-bright sources and the IRAC Peakers. 
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Fig. 4. Distribution of the 24/im sources selected as 
otically-faint IR-bright galaxies (blue dots) and with the 
BM/BX criterion (green dots) in the BzK diagram. Only 
sources at 1.5 < z < 2.8 are represented and galaxies sat- 
isfying the two selections are shown with pink color. The 
Optically-Faint IR-bright sources are characterized by red- 
der colors, most likely due to higher dust extinction. The 
red arrow represents the extinction vector. 



z ^ 2 (Figs. 2a & 3b) and at these redshifts the measure 
of 8/im luminosities (Ls/im) based on 24/im fluxes is there- 



fore almost independent of the SED templates assumed for 
deriving the /c-corrections. Also, the Sjam luminosity of star- 
forming galaxies mostly originates from the emission of the 
7.7/im and 8.6/im Polycyclic Aromatic Hydrocarbon (PAH) 
lines, which represent the most prominent features in the 
mid-I R spectrum of dusty galaxies ( e .g., Laurent et al. 
20001 [Dale et al.[ [2QQT| [Brandl et aTl |2006| [Smith et al. 



2007). These features are stochastically heated by the ra- 
diation field of young stellar populations, and as it was 
demonstrated by numerous studies of local sources, their 
luminosity is tightly correlated with the total Infrared lu- 
minosity of galaxies, and therefore with their SFR (e.g.. 



Roussel et al. 


2001 


Wu et al.| 20051 |Brandl et all |2006 


Calzetti et al. 


2007 


Diaz-Santos et al. , 2008 Goto et al. 



2011). Furthermore, observations of distant sources with 
Spitzer and more recently with Herschel have revealed that 
this correlation between the 8/im and the total IR luminos- 
i ty of star-forming galaxies extends up to at least z ^ 2 



Pope e t al.|[2008r|Rigby et al.||2008| 

^enendez-Delmestre et al. 1120091 E\h, 



^ ^_ [Bavouzet et al.||2008 

Menendez-Delmestre et al. 2009[ Elbaz et al., submitted) 
suggesting that the 8/im emission is a fairly good tracer of 
star-forming activity also in high-redshift galaxies. 

To convert the observed MIPS-24/im fluxes into lumi- 
nosities at 8/im we use d the library of IR SED templates of 
Chary fc Elbaz| ([2001). In this library the galaxy IR SEDs 
vary as a function of total IR luminosity (Ljr = i^8-iooo/xm) 
and for a given redshift the flux density measured at 24/im 
corresponds to a unique monochromatic luminosity at any 
IR wavelengths. However we stress again that for our cur- 
rent analysis the /c-corrections only depend on the shape of 
these templates between Aq = 8/im and Ai = 24/im/(l+z). 
At z ^ 2 they are thus barely sensitive to the choice of 
SEDs and do not depend on the uncertainties that have 
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been shown to affect t he extrapolations o f mid-IR fluxes to 



total IR luminosities (Elbaz et al. 2010). In fact we com- 



~i I I I r 



puted the rest-frame 8/im luminosities of the MIPS galaxy 
sa mple using the librari e s of IR SEP template s proposed 
by [Dale fc Helou] ( |2002D , [Lagache et aT] ( |2004D and |Rieke 
et aT ( |2009p , and we found virtuall y no difference with tfi e 
results obtained with the library of Chary & Elbaz ( 2001| . 
This effect will be discussed in more details in a forthcom- 
ing paper by Le Floc'h et al. (in prep.) on the evolution of 
the mid-IR luminosity function in COSMOS. 

Based on these Sjam luminosities we computed the lu- 
minosity functions (LF) associated with our different sub- 
samples over the 1.7 < z < 2.3 redshift range and we com- 
pared these results to the global LF obtained at the same 
redshifts from the fuh sample of MIPS sources in COSMOS. 
These luminosity functions are illustrated in the top panel 
of Fig [6} They were estimated using the 1/Vmax method 
(SchmmtJ 1968), which is advantageous in two ways: no 
fiypothesis on the shape of the LF is needed and the lu- 
minosity function is directly measured from the observa- 
tions. For each object we estimated the maximum comov- 
ing volume where it can be detected within our given red- 
shift bin (1.7 < z < 2.3) using Vmax — Vzmax ~ Vzmin^ 
where Zmin = 1.7 and z^ax is the smallest value between 
the redshift bin upper limit {z = 2.3) and the redshift until 
which the source could be observed considering the 24/im 
flux limit of our survey (i.e., F2A^m = 0.08 mJy). The un- 
certainties on the determination of the LFs were estimated 
from the Poisson noise associated to the number of sources 
detected in each bin of luminosity. Since we are only in- 
terested in quantifying the relative contribution of galaxies 
identified with the color selection techniques discussed ear- 
lier as a function of mid-IR luminosity, we neglected the 
effect of other uncertainties like the 24/im flux errors and 
the cosmic variance. 



5.2. Star Formation Rate Estimates and SFR density 

Because the mid-IR emission of star-forming galaxies orig- 
inates from dust features heated by ionizing photons com- 
ing from young stars, their mid-IR luminosity c orrelates 
pretty well with their star- formation rate (e.g., |Roussel| 



eraL][2QQTl|FQrster Schreiber et al.[|2QQ4l|Wu et al.l|2005p. 



For example Calzetti et al. " p007p and Dfaz- Santos et al.| 
( 2008 ) analyzed the mid-IR properties of HII regions taken 
from a sample of local star-forming sources. They estab- 
lished a significant correlation with a logarithmic slope ~ 1 
between their 8/im luminosity surface density corrected for 
stellar contribution and their luminosity measured from the 
Faa emission line corrected for extinction. Similarly, the 
combination of the mid-IR and far-IR observations of the 
nearby Universe conducted with IRAS, ISO, Spitzer and 
AKARI revealed a tight relationship betwe en Lsum the 



total IR luminosity of lo cal galaxies (e.g., Bavouzet et al. 



2008[ [Goto et al. 2011). The deepest observations of the 
sky undertaken with Herschel have recently shown that 
this correlation extends up to at least z ^ 3 (Elbaz et al. 
submitted) and it may thus represent a universal property 
shared among star-forming galaxies at all cosmic times. 

The relation obtained by Elbaz et al. (2011, submit- 
ted) between Lg^rn and the total Infrared luminosity is lin- 
ear {Liji/Ls^rn ^ 4.9 ± 0.2 dex). Assuming the standard 
conversion inferred by Kennicutt (1998) between the Lir 



COSMOS f24^„>80MJy with L7<z<2.3 



# Rodighiero +10 
■ Caput! +07 




ILO 1L2 1L4 1L6 1 L8 12,0 
Log(i/L (8Mm)) 




L5<z<2.8 
I ■ ■ ■ I 



1K0 11.2 



^^A 1K6 1K8 

Log(i/L (8Mm)) 



12.0 12.2 



Fig. 6. Top: The 8/im luminosity functions of MIPS galax- 
ies identified at 1.7 < z < 2.3 with the color selection 
techniques described in Sects. 3 & 4 (red: BzK; green: 
BM/BX; yellow: IRAC peakers; blue: Optically- Faint IR- 
bright galaxies) compared to the total 8/im luminosity func- 
tion measured in the COSMOS field (black solid line). 
Previous me asurements of the total 8/im LF at z ~ 2 



et al.j L2007j filled squares) are also shown for comparison. 



by [R odighiero et al. (2010, filled diamonds) and Caputi 

Error bars were estimated from the Poisson noise. Bottom: 
Fraction of galaxies identified with each of the color selec- 
tions as a function of Lg/im and in the 1.5 < z < 2.8 redshift 
range. Color coding is similar to that used in the top panel. 



and the star formation rate {SFR [M^yr ^] = 1.72 x 
10~^^ L I ji[Lq])^ we then obtain: 



SFRiMQyr-') 



8Atti X 10- 



(5) 



At the typical 8/im luminosities found in our high-redshift 
sample (i.e., Lg^m ^ 10^^ — lO^^L©), this is fully consistent 
with the relation we would get by combining the correlation 
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measured by [Goto et al. (2011) with the calibration from 6. Discussion 
Kennicutt||l998p : 



SFR{Me yr-') = (34 ± 9) x 10"^° x (Ls^^/Lo)^-^^ (6) 



Similarly, it agrees with the con version we wou l d obt ain 
by applying the relationship from Calzetti et al. (2007) to 
a typical star- forming disk (i.e., 2-5 kpc in diameter) and 
assuming a case B recombination for converting the Paa 
luminosity into a star for mation rate (SFR \Mc?)yr ~'^^ — 



6.79 X lO-'^^L paa \erg s-^] , lOsterbrock & Bochkarev 1989 



Alonso-Herrero et al.||2006p . This last approach relies on 
the assumption made for the size of the emitting region, 
but in the range of 8/im luminosities that we measured 
this effect is overwhelmed by the internal dispersion found 
in our sample and the uncertainty affecting the slope of 
the correlation between the Paa and the 8/im luminosity 



surface densities (Calzetti et al. 2007 Diaz-Santos et al. 



2008 see also Alonso-Herre ro et al.||2006[ ). This relation is 



also consistent with the Paa and the 8/im luminosities that 
were recen tly measured in a submillimeter lensed galaxy at 
z ~ 2.5 by Papovich et al. (2009). It suggests that it could 
still be applied to high-redshift sources even though it was 
derived from observations of the nearby Universe. 

Given the good agreement between these different rela- 
tionships, we thus inferred the total IR luminosity and star 
formation rate of the COSMOS MIPS sources using the 
linear relation obtained by Elbaz et al. (2011, submitted) 
as well as our Eq.[5) baring in mind that the other possi- 
ble methods would have led to very similar results. Based 
on these determinations of Ljr and SFR, we computed the 
contributions of the different populations of MlPS-selected 
sources discussed in Sect. 3 to the IR luminosity density 
of the Universe and the cosmic star formation rate den- 
sity in three redshift bins: 1.5 < z < 1.9, 1.9 < z < 2.3 
and 2.3 < z < 2.7. To reach this goal we derived lumi- 
nosity functions in a way similar to the method described 
in Sect. 5.1, and we integrated the LFs above the luminos- 
ity corresponding to the 24 /im flux limit of our survey 
(^24/im = 0.08 mJy) at the median of each redshift bin. 
Since we only aim at determining how the various galaxy 
sub-samples detected with MIPS contributed to the total 
SFR density, we did not assume any extrapolation of the 
MlPS-selected source population to the faint end of the 
luminosity function. Our results are shown in Fig.[7| which 
also illustrates the evolution of the total IR luminosity den- 
sity inferred by Rodighiero et al.| (|2010|) up to z ^ 2.5 and 
the cosm ic star formation history compiled by [Hopkins fc| 
Beacom ( 2006 J . The 1306 sources with no optical counter- 
parts are obviously not taken into account here but we al- 
ready saw that they contribute only a small fraction of the 
total number of galaxies identified in the COSMOS field. 
Furthermore, most of these unidentified sources have low 
24/im fluxes close to the sensitivity limit of the COSMOS 
MIPS data. This implies that their contribution to the IR 
luminosity density is even less than their contribution to 
the number of mid-IR selected sources, and therefore the 
lack of identification for these sources can not severely bias 
our results. 



6.1. Contributions to the IR luminosity function 

The top panel of Fig. [6] illustrates the luminosity functions 
associated with the MIPS sources identified with the four 
color selection techniques described in Sect. 3, compared to 
the total 8/im luminosity function observed at 1.7 < z < 2.3 
in the COSMOS field. The later is also co mpared with the 
8/im L Fs obtained by Caputi et al. (2007) and ,Rodighiero| 
etal] poTol ) from the GOODS and SWiRE surveys, which 
show relatively good agreement with our estimates. Except 
for the selection of Optically- Faint IR-bright galaxies, our 
analysis seems to indicate that the fraction of MIPS high- 
redshift galaxies selected with the BzK^ BM/BX and IRAC 
peaker criteria does not critically depend on the 8/im lumi- 
nosity itself. To better illustrate this result we reproduced 



the evolution of these fractions with L 



8/im 



in the bottom 



panel of the figure. To increase the statistics we extended 
the redshift bin up to 1.5 < z < 2.8, which corresponds 
to the largest redshift range where the efficiencies of our 
different color selections can be simultaneously compared 
with one another. We see again that the fractions of MIPS 
galaxies selected as BzK and BM/BX sources or as IRAC 
peakers are fairly constant over the range of 8/im luminosi- 
ties probed by the MIPS observations. They correspond to 
the average fractions we had already derived based on their 
number density and their redshift distribution. 

While the absence of correlation between Lg^^n and the 
fraction of BzK sources or IRAC peakers is not necessar- 
ily unexpected, the result that we find for the BM/BX 
galaxies may deserve further explanations. Indeed we at- 
tributed the bias observed in the BM/BX selection to the 
effect of dust extinction (see Fig. 2c). Given the global trend 
that exists between galaxy bol omet r ic lum inosities and the 
IR/UV luminosity ratio (e.g., [Belli [2003] ), one could have 
expected some broad correlation between the 8/im luminos- 
ity and the fraction of MIPS high-redshift sources missed 
by the BM/BX criterion. In the most luminous galaxies 
though, the reddening of the optical light does not corre- 
late anymore with the dust obscuration measured by the 



excess of Lir over Luv (e.g., Goldader et al. , 2002 Reddy 



et al., 2006). This can be understood, for instance, if on 
one hand the UV/optical emission originates from spatially- 
extended regions of star formation where dust grains fol- 
low a rather clumpy distribution, while on the other hand 
the IR light predominantly comes from very compact and 
optically-thick regions close the center of galaxies. Since 
the MIPS observations of COSMOS are only sensitive to 
the very bright end of the high-redshift galaxy luminosity 
function, we are likely in this luminosity regime where the 
E{B — V) extinction is not directly correlated with the level 
of star-forming activity, hence explaining the lack of trend 
between the IR luminosity and the detectability of galax- 
ies based on their UnGR colors. Finally, it is also possible 
that our result at the highest luminosities (Lg^rn > IO^^Lq) 
is slightly affected by some AGN contribution, which can 
become significant at bright 24/im fluxes (e.g., [Ma rtmez-] 
Sansigre et aI||2005l[Houck et al.| [20051 [Brand etal||2006| ). 



In particular, the blue continuum arising from unobscured 
type-1 quasars could potentially bias the BM/BX selection. 
However the majority of the most IR-luminous AGNs in 
COSMOS have also been detected in the X-rays ( [Hasinger 
et al. 2007 Brusa et al. , 2010). They have been removed 
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Fig. 7. Total Infrared Luminosity Density produced at respectively 1.5 < z < 1.9, 1.9 < z < 2.3 and 2.3 < z < 2.7 by 
24/im sources with F24/im > 0.08 mJy (black triangles), compared to the contribution of 24/im sources identified with 
the selection techniques described in Sects. 3 & 4 (red: BzK; green: BM/BX; yellow: IRAC peakers; blue: Optically- Faint 
IR-bright galaxies with photo-z identification). For galaxy samples where the MIPS 24/im emission is mostly produced 
by star formation (i.e., BzK and BM/BX soursces, IRAC peakers) t he IR Lum i nosity Density can be read as a Star 
Formation Density (right vertical axis) assum ing the conversion from i Kennicutt (1998). For comparison, the evolution 
of the total IR luminosity density inferred by [Rodighier o et al ^ ( |201Q| up to z 2.5 is illustrated by the light shaded 
region. The cosmic star formation history derived by Hopkins & Beacom (2006) is represented by the dark shaded area. 
The typical uncertainty affecting our estimates is shown by the verticaTTine in the top left corner of the diagram. 



from our initial sample and therefore they should not influ- 
ence the various trends that we find in Fig.|6] 

Contrary to the behavior characterizing the three rest- 
frame UV/optical selections, the fraction of Optically- Faint 
IR-bright galaxies is clearly rising with IR luminosity. This 
is consistent with what we had qualitatively inferred from 
the comparison of their redshift distributions at F24^m > 
0.08 mJy and F2A^m > 0.3 mJy (see Fig. 3). In the pre- 
vious section we argued that these sources are biased to- 
ward highly-obscured galaxies. Assuming that dust extinc- 
tion must correlate with the observed 24/im/optical flux 
ratio of galaxies at 2; ~ 2, it could thus explain the trend 
that we find with luminosity. Furthermore, this trend could 
partly originate from an increasing contribution of obscured 
AGNs to the mid-IR lu minosit ies of galaxies as a function 
of 24/im flux (e.g.. Brand et al. 2006^ Me nendez-Delmestre] 



et ah] 2009). This rising contribution would boost their 
mid-IR emission with respect to their optical luminosity, 
hence leading to an r+ — [24] color satisfying the selection 
criterion of this population. 



6.2. Implication on the cosmic star formation history 

Even though we did not account for the contribution of 
galaxies below our 24/im flux limit we see on Fig. 7 that 
star-forming galaxies detected in the COSMOS mid-IR sur- 



vey contribute a substantial fraction of the total IR lumi- 
nosity de nsity at z ^ 2, in agreement wi th previous fi nd- 



ings (e.g., Caputi et al. 2007; Rodighiero et al. 2010). It 



implies that the incompleteness which affects the different 
color selection techniques previously described can result 
in non-negligible biases with respect to the global picture 
of galaxy evolution inferred from complete samples of high- 
redshift sources. For example we find that the MIPS sources 
selected with the BM/BX criterion contribute only ~50% 
to the IR luminosity density produced by the galaxies with 
> 0.08 mJy at z ~ 2. Similar to what is observed for 
the other color selections, this fraction is roughly constant 
over the whole 1.5 < z < 2.8 redshift range. In the case of 
the BM/BX technique we attributed this bias to the effect 
of dust extinction, and since the obscuration globally de- 
creases with bolometric luminosity (Reddy et al., 2006 J the 
global incompleteness affecting this selection could proba- 
bly be smaller if we had considered a deeper 24/im selec- 
tion. However, it is interesting to note that a very similar 
result was also independently obtained by |Guo fc White] 
(2009) based on a semi-analytical approach and using the 



updated ver sions of the Millenium Simulations (Springel 
"2005|). According to the models they found that a 



et al. 



large fraction of the total star formation density at z ~ 2 
originates from galaxies that are too heavily obscured by 
dust to be selected with the BM/BX criterion. In princi- 
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pie, incompleteness corrections can be applied to account 
for the difference between the properties of galaxies satis- 
fying this color selection and the properties characterizing 
the global underlying population of star- f orming sources 



at these redshifts (e.g., Reddy et al.[ 2QQ8t. However, our 



analysis and the results from Guo & White (2009) reveal 



that a substantial extrapolation accounting for half of the 
most actively star-forming sources (i.e., Lg/^m > lO^^"^^©^ 
corresponding to SFR > lOOM0yr~^) would be required 
to also take into account the large contribution of dusty 
luminous galaxies to the cosmic star formation density at 
1.5 < ^ < 2.5. 

Similarly, we find that if the IRAC peaker selection is 
probably well suited to the identification of galaxies char- 
acterized by a pronounced stellar bump, it can clearly not 
provide a complete view on the total contribution of mas- 
sive star-forming sources to the build-up of stellar mass 
at high redshift. Over the whole redshift range where this 
selection technique can be applied (1.5 < z < 3) up to 
~30% of the total 8/im luminosity density produced at 
the bright end of the luminosity function is contributed 
by massive sources escaping the IRAC color criterion previ- 
ously discussed. In the case of the Optically- Faint IR-bright 
sources we find that this fraction can even go up to ~70%. 
Although the incompleteness affecting these two selections 
could be partly explained by the hot dust continuum of 
dust-obscured AGNs contributing to the rest-frame near- 
IR emission of the most luminous (i.e., I/gum > lO^^Lc? 



galaxies at high redshift (e.g., |Houck et al.[ |2005[ [Desai 
et al. , 2009), it is unlikely to be the case for the majority 
of distant sources with more typical luminosities and which 
are known to exhibit strong PAH features ch aracteristic of 



star-f orming activity at mid-IR wavelengths ( Fadda et al 
2010). Furthermore we recall the strong correlation that has 



been recently established between the 8/im luminosity and 
the total IR luminosity of galaxies up to 2; ~ 3 (Elbaz et 
al., submitted), which reinforces the idea that the mid-IR 
luminosity characterizing the general galaxy population is 
mostly powered by star-forming activity. We thus conclude 
that a non negligible fraction of the cosmic star formation 
density is missed by these two selection techniques. 

On the other hand we note that the BzK approach en- 
ables an almost complete selection of star-forming galaxies 
at z ~ 2 at least in the range of luminosities that can be 
probed by our MIPS observations. This high level of com- 
pleteness is observed not only for the bulk of the sample but 
also for the most luminous - and potentially the most highly 
obscured - galaxies. Using the optically-selected catalog of 
photometric redshifts in COSMOS ( Ilbert et al. , 2009) we 
verified that the reliability of the BzK color criterion to 
select star- forming galaxies in the 1.4 < z < 2.5 redshift 
range remains also valid at fainter luminosities. It suggests 
that this technique provides a quite powerful way to identify 
galaxies at this epoch of cosmic history, independently of 
the possibility to determine accurate photometric redshifts 
for these sources. 

More globally speaking, the fact that a non negligible 
fraction of high-z sources can be missed by optical color 
selections has also been observed by other groups using dif- 
ferent kinds of analysis, and this bias likely originates not 
only from dust extinction but also from the broader diver- 
sity of high-redshift galaxy properties than was assumed 
in the definition of the color criteria that were explored so 
far. For instance, Ly et al. (private communication & in 



prep.) show that surveys of [OH] emitters can pick up a 
non negligible fraction of z ~ 1.5 galaxies that are missed 
by th e BzK and BM/BX techniques. Similarly, Caputi] 
et al.| (|2011| ) found that at z = 3 - 4, -30% of the mas- 
sive galaxies (M>10^^ ^0) would be missed by deep optical 
surveys, underlying the importance of deep near-IR surveys 
to supplement the identification of high-z galaxies affected 
by extinction. Based on deep s pectroscopy taken as p art 
of the Vimos VLT Deep Survey, |Le Fevre et°aL| ( |2005[ ) fi- 
nally underlined that magnitude-selected samples provide 
a much more complete census of the high-redshift galaxy 
population than color selection techniques. All these results 
converge in the way that the identifications of distant star- 
forming galaxies based on UV/optical color criteria result 
in a view on high redshift star formation that can not be 
fully representative of the global picture of galaxy evolu- 
tion. 



7. Summary and conclusions 

Using the exquisite multi-wavelength imaging obtained 
in the COSMOS field we analyzed the broad-band 
optical/near-IR colors and the photometric redshift distri- 
bution of a complete sample of distant (1.5 < z < 3) lumi- 
nous galaxies selected at 24/im with F24^m > 0.08 mJy. 
To quantify how selection effects and incompleteness is- 
sues can bias the different color selection techniques of- 
ten used for identifying populations of high-redshift star- 
forming sources (e.g., BzK, BM/BX and "IRAC peaker" 
selections, "Optically- Faint IR-bright" objects), we mea- 
sured the fraction of MIPS-24/im galaxies satisfying these 
color criteria within the redshift range where each of them 
is expected to apply with high efficiency. We also derived 
reliable estimates of their rest-frame 8/im luminosities with 
minimal effects from the /c-correction, which allowed us to 
constrain their star-formation rate and their contribution 
to cosmic star formation density. Our results can be sum- 
marized as follows: 

1. The BzK color criterion allows an almost complete 
(~90%) identification of 24/im sources at 1.4 < 2; < 2.5. 
The rate of these identifications does not change with 
IR luminosity, which suggests that this technique 
provides a highly reliable way to identify distant 
galaxies independently of obscuration. This can be 
explained by the weak dependence of the BzK variable 
{BzK = [z - K] - [B - z]) on the effect of dust 
extinction at 2: ~ 2. However we emphasize that the 
success rate of the BzK selection strongly relies on 
the availability of deep imaging at optical and near-IR 
wavelengths. With a cut at Kg = 22 (AB), the contri- 
bution of Bz K-identifi.ed galaxies to the IR luminosity 
density produced by sources with F24^rn > 0.08 mJy 
would fall down to ~40%. 

2. We adapted the original BM/BX criterion to the 
existing COSMOS optical photometry using the 
deep imaging obtained with the ii*-band, Vj-band 
and z+-band filters. We found that up to 50% of 
the MIPS selected galaxies at 1.5 < z < 2.8 are 
missed by this selection technique. We explain this 
result by the effect of dust extinction, which more 
strongly affects the rest-frame UV emission of galaxies 
and makes their SED redder than the typical blue 
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colors characterizing galaxies satisfying the BM/BX 
criterion. T his is co nsistent with the results recently 
obtained by Guo fc White (2009) b ased on Millenium 
simulations ( [Springel et al!] ~2005), who found that 
a large fraction of the cosmic star formation density 
probably originates from massive sources that are 
too heavily obscured by dust to be identified with 
the UV color selections. At the bright end of the 
bolometric luminosity function, the extrapolations and 
incompleteness corrections usually applied to account 
for this bias must be therefore quite large and uncertain. 

3. Up to ~30% of the IR luminosity density produced at 
z ~ 2 by sources brighter than 80/iJy originates from 
galaxies that fail to be identified thanks to the spectral 
shape of their stellar bump feature at rest-frame 
1.6/im (i.e., the so-called ^^IRAC peaker^^ selection). 
For the brightest galaxies in the mid-IR, this effect 
could be explained by an important contribution of 
dust-obscured AGNs to the near-IR emission observed 
in the IRAC bands. For the bulk of the population 
though, it reflects the larger uncertainties affecting 
the IRAC photometry as well as the wide diversity of 
near-IR SEDs characterizing dusty galaxies, due to the 
combined contribution of hot dust and stellar emission. 

4. The selection of Optically- Faint IR-bright sources based 
on the R — [24] color is meant to identify dusty galax- 
ies at 1.5 < z < 3. While this technique is efficient 
at the very bright end of the mid-IR luminosity func- 
tion where galaxies are characterized by extremely high 
24//m/optical flux ratios, the fraction of sources selected 
with this criterion rapidly decreases at fainter luminosi- 
ties. Assuming the depth of the COSMOS MIPS survey 
{F24f,m > 0.08 mJy) we find that only -25% of the IR 
luminosity density of the Universe at z — 2 is produced 
by galaxies with such extreme colors. 

These results suggest that color selections of distant 
star-forming galaxies can be affected by substantial biases 
and incompleteness, and therefore they must be used 
with strong caution. Depending on the spectral features 
or the SED range that they are supposed to probe, 
these high-redshift selections are probably less subject 
to contaminants compared to the use of photometric 
redshifts with modest quality if the later are strongly 
affected by catastrophic failures at z > 1. However we 
have shown that the identification of distant galaxies 
solely based on optical/IR color criteria will most often 
provide an incomplete view on the whole population of 
high-redshift star- forming sources, hence requiring large 
and uncertain extrapolations to account for their incom- 
pleteness. Complete identifications based on high-quality 
spectroscopic redshifts or accurate photometric redshifts 
thus appear the most reliable approach for probing the 
formation and the evolution of galaxies in a global and 
cosmological context. 
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